1 



SINP-TNP/06-23 

Quantization and Conformal Properties of a 
Generalized Calogero Model 

S. Meljanac ^, A. Samsarov ^ 
Rudjer Boskovic Institute, Bijenicka c.54, HR-10002 Zagreb, Croatia 



B. Basu-Mallick^, Kumar S. Gupta ^ 
Theory Division, Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Calcutta 700064, India 



Abstract 



We analyze a generalization of the quantum Calogero model with the underlying conformal symmetry, paying special 
attention to the two-body model deformation. Owing to the underlying SU{1, 1) symmetry, we find that the analytic 
solutions of this model can be described within the scope of the Bargmann representation analysis and we investigate 
its dynamical structure by constructing the corresponding Fock space realization. The analysis from the standpoint of 
supersymmetric quantum mechanics (SUSYQM), when applied to this problem, reveals that the model is also shape 
invariant. For a certain range of the system parameters, the two-body generalization of the Calogero model is shown 
to admit a one-parameter family of self-adjoint extensions, leading to inequivalent quantizations of the system. 

PACS number(s): 02.30.Ik, 03.65.Fd, 03.65.-w 

Keywords: Bargmann representation, ladder operators, shape invariance, self-adjoint extensions 



^ e-mail: meljanac@irb.hr 

^ e-mail: asamsarov@irb.hr 

^ email: biru@theory.saha.ernet.in 

* emaihkumars.guptaOsaha.ac.in 



2 



INTRODUCTION 



The structure and apphcation of the Calogero model 0, Q and its various descendants is a subject that receives 
much attention. Since these are examples of man y-body ex actly solvable models which appear in various contexts in 



physics as well as in mathematics |,'Ul4ll5ll(ll7ll8ll9l llflllllll2lll,l |. it is of considerable interest to find its generalizations 
which are exactly solvable and integrable. In particular, it is appealing to investigate whether some sort of modification 
of an exactly solvable model will affect its integrability. The dealing with such many-body problems in one dimension 
has been shown up as particularly advantageous since there exist several algebraic techniques that are applicable in 
this case, due to highly restrictive spatial degrees of freedom. It is also interesting to find out all possible boundary 
conditions that render the Hamiltonian of the system self-adjoint [l^ and to analyze the nature of the corresponding 
spectrum. 

In this paper we investigate one special class of deformation of the quantum Calogero model. As a prototype, we 
study the two-body model as a particularly convenient one for a complete elaboration of the dynamical structure of 
the problem, which is basically the same as for all many-body problems with the underlying SU (1,1) symmetry. This 
specific case will be approached from more different directions that include the Bargmann representation analysis, the 
ladder operators formalism and the SUSYQM analysis which is based on the shape-invariance property of the model 
in question. We also study the self-adjoint extensions of the two-body model and show that for certain values 
of the system parameters the model admits a one-parameter family of inequivalent quantizations. At the end, the 
analysis will be expanded to include the A^-body case as well, with the result for the complete spectrum of the 3-body 
case stated explicitl y. T he emphasis will particularly be given to the Bargmann representation approach which was 
put forward in Ref. [l5j in context with the general method for integration of the multi-species and multi-dimensional 
generalizations of the Calogero model. This approach relies heavily on the conformal invariance of the model under 
consideration. 



2. BARGMANN REPRESENTATION ANALYSIS 

In a recent paper jl5j . a general procedure for integrating a many-body quantum systems with the underlying 
SU{\, 1) symmetry was set up. This procedure is based on the fact that the N— body quantum system, possessing 
conformal symmetry, can be mapped onto the set of N harmonic oscillators in arbitrary dimensions and with a 
common frequency lo. In this paper we investigate the system described by the Hamiltonian 

„ 1 / 92 92 \ 1 2, 2 2^ ^ ^Ji 

as a simple example of the model and techniques that were put forward in The above Hamiltonian clearly 

represents two interacting particles in one-dimensional space bounded by a harmonic force. For simplicity, masses of 
the particles are set equal to 1 and ?i = 1. It is convenient to introduce the following set of operators: 
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7"- — 2 'n,i=l ^ V{xi,X2) — 2 (ar2 + r Or + ) V, (2) 



2 2 Or ' 2' 



2{xi-X2f ' 2(^2 +2:2) ■ 



v = - ^ + :T7::3^- (3) 



with V given as 



If we set iJ, = —\ in , we arrive at the model considered in Ref . [ig. The operators T+,T-,To in are expressed 
in terms of the polar coordinates xi — r sin (j), X2 — r cos 4>. A brief inspection of the Hamiltonian |^ reveals that it 
belongs to the class of conformally invariant systems ,15j , with the potential V being a real homogeneous function 
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of order —2, i.e., satisfying the relation 



i—l 

Indeed, it is straightforward to show that the generators Q satisfy the SU(1,1) conformal algebra 

[T_ , T+] = 2To, [To, T±] = ±T± (5) 

and that the Hamiltonian (QJ can be represented in terms of them as H = — T_ + oj^T^. The problem of solving 
for the eigenstates and for the spectrum of the Hamiltonian Hijjn^k = En,ki^n,k, will now be transferred into 
the two-oscillator eigenvalue problem for the Hamiltonian 20-" To, which is usually referred to as a transition to a 
Bargmann representation. This can be achieved by applying the transformation 

H = 2ujSToS-\ (6) 

where 

S^e-"^+e-i=^~. (7) 

After carrying out the transition to the Bargmann representation, we have to solve the eigenvalue problem for Tq, 
with an additional constraint 

r_A„ = 0, roA„ = ^A„, n>0. (8) 

As it stands, there are many eigenfunctions of Tq. In fact, every homogeneous function, no matter whether it is 
rational or irrational, is an eigenfunction of Tq. However, among all of them we have to pick up only those that are 
annihilated by T_ . These form an infinite number of vacua A„ , upon which the equidistant towers of states are 
built up, with an elementary energy step 2uj between any two neighbouring states. The ground state -00,0 of Q 
is required to be a square-integrable function. This will be the case if the ground-state e nerg y ujeo, Cq > 0, is 
higher than ^, the condition which is connected with the existence of the critical point [tR lisf . Then we can write 
"00,0 = 'S'Ao, where Aq is a homogeneous function of the lowest degree and of the lowest energy weo- The other vacua 
A„, n > 0, of higher degrees of homogeneity and with energies we„, are also mapped to ipn,o = SAn- As far as 
the relations (O are concerned, the generators ((SJ can be viewed as creation and annihilation operators acting on the 
eigenstates of To. In respect of this, the excited states Vn./c of ^ can be constructed as 

^nM^ST+'^An, fc = 0,l,2,...; n>0, (9) 

with energies 2w(fc -I- and S given by Q. Thus, to conclude, we have all solutions grouped into equidistant 
towers of states based on SAn, n > 0, and for a given n > 0, the spectrum is equidistant with an elementary step 
2w. 

It may be noted that the Casimir of the SU{1, 1) algebra © is given by 

C =T^ -To-T+T-, (10) 

which commutes with the Hamiltonian, the latter being an element of the algebra (O. The given system is thus 
integrable with the Hamiltonian and the Casimir as the two conserved quantities. By using ((HJ, © and l(TT)| . we see 

that 

Thus it is evident that each tower of states built on the vacuum A„ provides an irreducible representation of the 
SU{1, 1) algebra classified by the eigenvalue of the Casimir C given in 
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The aforementioned procedure will now be applied to the model Hamiltonian |^ to yield the eigenstates and the 
spectrum. First, we have to solve the equation T_A„ = 0, n > Q, for the potential 
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y2dxl 2dxl 2{xi-X2f 2(.x2+x2)^ 
In polar coordinates this equation reads 



A„=0. (12) 



- 4 A„ = 0. (13) 



Qj,2 J, Qj, J.2 Qfjp. r2(l — sin20) r2 
If we follow the lines of Rcf. , by separating the variables as A„ = M(r)$(0), (I13|) reduces to the pair of equations 

d'^u 1 du C + 



.92$ / A 



-u = 0, (14) 



2 V 1 - sin 20 



-C $ = 0. (15) 



The first one integrates to u ~ j.^/c+t'.^ while the second one appears to have physically acceptable solutions only 
for the special values of the coupling A and the separation constant C. In order to find these solutions, we skip to 
the angular variable = <^ — ^, running through the interval Q < < t:, and then factorize the function $, 

$„(^) = (sin0)''/„(cos0) = {^xneyUx). (16) 

where the variable x = cos 9 has been introduced and n is the quantum number labeling the vacua A„ . After 
inserting (|16|l into equation H15() . we get the equation for the functions /„, n > 



, , sin^6i - {2v + l)x^!-^^swL^e + (v{v - 1) - -)x^ - (v + - - C)s\n^e fJx) = 0. (17) 
ax'' ax \ 2 2/ 

Physically acceptable solutions to the above equation emerge for the special choices of the coupling and the separation 
constants A and C, respectively, 

A = 2i^(i^-1), C^{n + iyf, (18) 



in which case Eq. (|17|l becomes the equation for the Jacobi polynomials Pn""'''^ 

(1 _ ^2^f^l^ _ (2^ + l)a;%M + n{n + 2v)Ux) = 0. (19) 



Thus, as a solution to itli^ we have /„(x) = Pn {x), with a — ^ . Altogether, this gives us the following 
expression for the vacua in the Bargmann representation: 



A„ = ^ r\/^(sin(?)7„(x) = rV(^^+^^(sin^)"P,('^''^)(cos^), a = (20) 
In view of eqs. (jSJ, the excited states in the Bargmann representation are obtained as 

T+''^A„ = (-r2) rV(«+'^)'+f(sin0)"F,('^'°)(cos0). (21) 

With the help of the transformation 10, we now transfer the results back (see eq. Q) to the original problem to 
obtain the eigenstates of Q, 



V-n.fe = ST+'^A,, ~ e-^+ -ir^T^-^ (^^V(«+-)^+f +2^ (sin 0)'^P,(-'-) (cos 0)j . (22) 

1=0 ' 
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For the moment, we shall not border with the proportionality constants because they are irrelevant for the discussion 
below. This is the reason why the sign of similarity appears in (|22|l and in the majority of subsequent relations. Later 
on, when we consider a construction of the bosonic ladder operators, we shall have to take care of the normalization 
of the wave functions and every single proportionality factor will be important. 

In calculating ( I22|l . we shall make use of the fact that the successive application of the operator T_ to the states 
(|21|l leads to an expression of the form 

tJ (t+'^A,,^ - 2"'2fc (2k + 2^J{n + uf +f?j {2k - 2) ^2fc - 2 + 2\J{n + vf + i^j x 
{2k - 4) ^2/s - 4 + 2\l{n + vf + 1^ ■ ■ ■ {2k ~ 2{l - 1)) ^2fc - 2{l - 1) + 2\l{n + vf +1^ x 



rVM^^^)^+2fe-2i (gj^ ofpM (cos e). (23) 

When I is equal to fc, the above expression is, essentially, the Bargmann vacuum A„, up to some constant. Due to 
the first of Eqs. the next application of r_ yields zero, causing the series in H22|l to terminate, T_' ^T+'^An^ = 
for I > k. Relation f 122(1 now becomes 

~ 2-'=e-^'^VV("+'^''+'^(sin efP^;'^''^ (cos 9) x 
ofe+l / I \ ofc+ 



2^2'=-^— fc [k+J{n + vf + /i r^'^-^^ ^k{k-\) k + V ^vf ^ A[k~\^ J {71 + vf + ^i] r^""-^^ 



2uj 



(2w)"2! 



+ (-!)'' ■^^^^fc(fc - l)(fc -2)---2-\[k+ Y^(n + vf + (^fc - 1 + + j/)' + • • • (^1 + Y^(n + + M 



= ,-.,-t.V^/(^^^^(sin0)''p(-'')(cos 0) ^ (-1)'^^)'=-' Q ^1 + V^^^^ , (24) 



(=0 

where the binomial coefficients in (|24|l have the following meaning: 

a\ _ r(a + l) 



/3!(a-/3)! r(/3 + l)r(a - /3 + 1) ' 

with r being the Euler gamma-function. 

If we rearrange the sum in (|24|l . we finally obtain the result 



[k + ^^{n^vf ^ ! 



= (-l)'=fc!c^-^e-t-VV(«+-)^+M(sin0)-p/^-.-)(cos0)L/^"+"^'+^(L^r2) (25) 

for the eigenfunctions of the Hamiltonian (Q. In obtaining this result, the sum in the first line of equation H25|) is 
recognized up to the irrelevant factor (anyway, it enters the overall normalization constant of the wave function), 
as the power-series expansion for the associated Laguerre polynomials. 



^ (_1)™ /fc + a 



m=0 
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Considering the spectrum, as it is unaltered by the transformation (j?)), we conclude that the eigenenergies of Q are 
nothing else but the degrees of homogeneity of the states in the Bargmann representation. When To is applied to the 

Bargmann vacua fiOfl . we obtain ^ {n + v)^ + ^ + 1^ A„ and according to the second one of equations (|SJ|, this 

has to be equal to ie„. Knowing the dynamical structure of the problem considered ( see Fig. 1), and anticipating 
that for a particular vacuum labeled by n, the energies of the neighbouring excited states are separated by 2uj, the 
spectrum is easily found to be 



En,k = (^{en + 2k)=uj (^{n + vf +11 + 1 + 2fc j , 



(27) 



with V being determined by (|18|l . 



3. LADDER OPERATORS 



As was said previously, in order to find bosonic ladder operators related to the eigenfunctions we were 

in search for, the knowledge on the detailed form of the normalization factor is essential. Thus, after we utilize 
the normalization properties of the Jacobi and Laguerre polynomials, we are left with the normalized version of the 
eigenfunctions (|25|) of the model Hamiltonian J^l, 



\ 



-^ ^^l^^^'kl n^[n + .)rin + 2.) ,-f.VV(^^(3in.)V(-)(cos^)Lf^K^) 

(28) 

For brevity, we shall omit the tilde-symbol from the wave functions and in the whole subsequent exposition it will be 
understood that they are normalized to unity. As we are searching for the operators that provide transitions between 
ground-state configurations of neighbouring towers in Fig. 2, we basically look for the operators connecting Jacobi 
polynomials of successive order. A straightforward way of deducing on the form of these ladder operators (in Fig. 2 
they are designated by B and B^) is to invoke the recursive relations 19] for Jacobi polynomials, 

(2n + a + b){\ - x^)^Pl^^''\x) = 



n[{a - 6) - (2n + a + 6)2;]/^^'^'^) (x) + 2{n + a){n + 6)1^^^'^^ (x), 



(29) 



2{n + l)(n + a + h+ l)(2n + a + h)P';;X^^ {x) = 

{2n + a + b+ l)[{2n + a + b){2n + a + b + 2)x + - b^]P';°-^\x) - 2{n + a){n + b){2n + a + b + 2)P^"'J^(x). (30) 

Using the recursive relations H29|l . (|3U|1 . we find recursive relations for the normalized ground-state energy eigen- 
functions V'n.Oj 

[{n + v)x X^)-^]i>nfi = 

ax 



{n + a) 



J{n + i^y+fi + l r(n + j/+ 1/2) 



(31) 



[{n-iy + 2a+ l)x - (1 - x^)-j^]i^n,o 
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(n+ l)(n + 2a+ 1) 
n + a+1 



\ 



4— ir [J{n + vY + n + l]T{n + v + 1/2) 



-e 2 ■ 7- 



(i-x^)-^pr+r(x), (32) 



where x = cos6 and a = — 5. 



n = 



n= 1 



n = 2 



n = 3 




Figure 1: Horizontal shifts in the Fock space of states are accomplished by the operators -B^ and B, whereas the 
vertical shifts are provided by the pair of operators A2 and A2 ■ 



The ladder operators b and 6^ that shift the neighbouring vacua into each other can be read out from the above 
recursions 



6=[i(JV + >/) + (l-i^)— ] 



i 



^^(N+u-if+t, {N + v-l) 



r 



V(w+i')'+M {N + v){N + 2v -I) 



r (^^J{N + v-if + iJL + i^ 



(33) 
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= [x{N + x"^) 



dx 



with N being the number operator defined as 



N 



ii-xr^'pt-\x) 



\i-xr^'pt^\x) 



The exphcit form of this operator can be found with the help of Eq. p5l ) and looks like 

, A 



N = VL2 - 



L2 = - 



1 - cos 20 ■ 

Since N is only angular dependent, relation H35|) can be straightforwardly extended to 

Straightforward calculation shows that the ladder operators b and are bosonic, 

[b,b^] = l, 

together with 

[N,b]^~b, [N,b'i]^b\ 



(34) 

(35) 

(36) 
(37) 
(38) 
(39) 



resulting in the simple relation including the number operator N = Wb. However, they are still not the ladder 
operators we are looking for because they do not shift the neighbouring eigenstates between each other, as it is readily 
seen from (|31|l and l|32|) . To overcome this problem, it is convenient to consider a certain type of operators that are 
realized by means of the similarity transformation applied to the bosonic operators b and 6^ 



5t = rx/I^5t^-x/^^_ 
These operators have the desired properties, namely, 

i.e., they are bosonic, satisfying 



(40) 
(41) 

(42) 

(43) 



Owing to the similarity-transformation kind of relation between the operators 6,6^ and B,B^, the latter pair 
retains the simple relation to the number operator, namely, N = B. Having established the form of the operators 
(|40|1 and (|41|l . we have succeeded to describe the horizontal shift in the Fock space of states, depicted in Fig. 1 by 
horizontal arrows. The vertical shift in Fig. 1 still remains to be described. However, it is easily accomplished by the 
transformed Bargmann-representation creation and annihilation operators |(5J), 

Af = ST±S-' = 1-(ujT+ + -t\ t To, 



(44) 

^ \ ^ / 

where T+,T_,To are conformal generators |(2J) and S is the transformation Q. Of course, the state of the lowest 
energy |0) = ipn=o,k=o is annihilated by both of the operators B and A^. Now, the general Fock-space state 
of Fig. 1 is easily obtained by the successive application of the ladder operators and to the vacuum state 

|0) = ■>pn=0,k=0 



= (A+)'(i?t)"|o) 

and in the coordinate representation it is explicitly described by the expression H28|l . 



(45) 
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4. SUSYQM ANALYSIS OF THE RADIAL PART 

The model under consideration obviously possesses some supersymmetric features and can be approached from the 
point of view of the supersymmetric quantum mechanics |2,'i 124. l25l | . To see this, we start with the Hamiltonian Q 



By making a transition to polar coordinates and by separating the variables in the manner as 0) = u(r)<i>(0), 
the Schrodinger equation H'ii = E'^ becomes 

I ^d^U I du 2 4 nr. 2 1 9^$ ^ 

u or^ u or 1 — sm 20 <3> c/ip-^ 



where C is the separation constant. This separates to the following two equations: 

ld_ 2 2 

Sr^ r dr r 

A 



92 Id 2 ? C + /i\ , , 

Zir+^'r' + —JT^]^-Eu, (48) 



9(/)2 l-sin2</> 
By introducing a new function V-'(^) = V^u{r), ()48|l simplifies to 



$ = C$. (49) 



1 _ 22 + - 3 



2V a,2+^^ + — f^)V^ = ^V'. (50) 



From this equation we can immediately read out the radial Hamilton 

where the parameter a = \JC + /i is introduced, with C determined by H18|l . The radial Hamiltonian H51|l is shape 
invariant and this property can be verified by introducing the corresponding superpotcntial 

a + ^ 

U^ujr ^, 52 

r 

and by factorizing the Hamiltonian (|51|) with the help of the operators v4(°\ A^'^^'^ introduced in the following way: 



_ , ^ . - , - ^ I J ' ^) . (54) 

dr ) dr r ) 

Now we have 

i/,, = A(")^A(")+eo, (55) 



where eg is given by (|27|l as cq = £'n,o = ij-'(q: + 1). Let us now create a supersymmetric partner Hamiltonian h'^^ 

ijW =A(o)A(")Veo. (56) 



of the Hamiltonian = Hf^ , 
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'Pn,k 



'f^n,fc-l 



(0) 



(0) 

,fe-i 



(0) 
n,2 



(0) 
rt,0 



-^(0)1- 



^^^^ 



(1) 

n.fc — 1 



fe-2 



(1) 




n.O 



(2) 
n,0 



n,0 




77. = 



ij(i) 



7J(2) 



Figure 2: A pattern for constructing the eigenstates of the model Hamiltonian H from the eigenstates of the 
supersymmetric partner Hamihonians H^'^K 



The factorization of ( I56|l in a way pursued in H55|l can be achieved by introducing the operators 



3 



^/2\dr r j 

It is easy to check that the two sets of operators, and (|^. 

condition: 

^(o)^(o)t^^(i)t^(i)^g^^ (59) 
where ei = 2a;. 

A somewhat more detailed insight into the relationships among the quantities we are dealing with here is obtained 
by taking a look at Figs. 1 and 2. From Fig. 1, let us take an arbitrary tower of states, which, for the sake of 
argument, we can take to be the tower built up upon the vacuum state labeled by n. This tower coincides identically 
with the zeroth tower shown in Fig. 2, with other towers in Fig. 2 being the eigenstates of the corresponding 
supersymmetric partner Hamiltonians of the Hamiltonian (|55|l . This correspondence is made obvious (see Fig. 2) by 
labeling each radial wave function by superscripts in parenthesis. In this way, the superscript i put on a particular 
radial function designates that this function is the eigenfunction of the super partner Hamiltonian H). , and so on. 



(58) 

satisfy the SUSYQM shape-invariance 
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Each supersymmetric partner Hamiltonian Hr"^ has its own set of eigenstates ipn\ spectrum eI^\, 

with n, i fixed and k being a nonnegative integer, 

^^^VSl = <ltfl k = 0,l,... (60) 

Note that the states from different towers but in the same horizontal line have the same energy and are transformed 
into each other by a simple pattern 01 • 

Having the whole picture set up, the basic vacuum ( the state V'i°o) obtained from the condition that the operator 
A'^^') should annihilate it, 

A^'^^Z^o - 0. (61) 

Equation integrates to 

,pl^l^r''+h-^^\ (62) 
In a similar way, the vacuum state of the i = 1 tower is obtained as 

^^'VSLo = 0, (63) 

leading to 

V'2Lo--"+'e-^'■^ (64) 

The first excited state V'i''l;=i and the corresponding energy E^l,^^ of the original Hamiltonian (|?T|l follow as 

<Li-^^"^Vi:Lo, (65) 

giving 

V-i^Li -e-t'-V^+^ft^r^- (a + 1)] =e-^''V"+5L5'(wr2), (66) 

<Li=<Lo = eo + ei=eo + 2c^, (67) 

where L^(a;r^) is the associated Laguerre polynomial and the last equation follows from the combined application 
of ISU, and (EOl- 

By following the same line, we are led (see Fig. 2) to a pattern for constructing a general radial excitation together 
with the corresponding excitation energy. First, we introduce the pair of operators 

^I.^t I f d a + k+^\ 

The conjecture, now, is that after applying H69(l j times successively, 

we get 



s=0 
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where the variable z = lur^ has been introduced. This expression can be proved by induction. For j = 1, the 
expression (|7(J|I reduces to 

which coincides identically with V'i'T^^ ~ ^^*'~^''^V'S- ^or the step of the induction, let us assume that <(7n|l holds 
for J. Then we have 



^ , , , wr ^ r ' 2 e 2 ■ r 

V2\dr 



By rearranging the sums and by using the properties of the binomial coefficients, the expression written above becomes 



J / t~; \ s / \ s 



s=0 



and this is exactly H70|) for j ^ j + Now that we have verified relation H70() . the eigenstates of the radial Hamiltonian 
= Hr^'' follow directly from (|70|l by setting j = k, 



s=0 



^fc\ , f a + k 



By rearranging the sum in the above expression, (|71|l is recognized as an expansion (|26f) for the associated Laguerre 
polynomials, 

<i~r-"+^e-^'^^L^K2), (72) 

as it should be, to coincide with the results obtained previously by pursuing other methods. To complete the SUSYQM 
analysis, we turn onto the spectrum of the model. This analysis, of course, has to yield the same result as obtained 
before. While carrying out the consideration, one should note that the shape-invariance condition (|59|l holds at the 
general level, i.e., 

^(fe-i)^(fc-i)t ^^(fc)t^(fc) ^^^^^ efc = 2w, fc = l,2,3,..., (73) 
as well as that the operators A^*^^ annihilate the corresponding vacua 

AWvS = 0, fc = l,2,3,... (74) 
Also, as Fig. 2 suggests, the energies of the states in different towers, but in the same line are the same 

p(0) _ p(l) _ p(2) _ _ p(fc) 
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Following the construction pattern for the supersymmetric partner Hamiltonians Hr , their form is immediately 
deduced from the shape-invariance condition (|73|l 

fc 

=^(fe)^^W^^g^.^ fc= 1,2,3,... (76) 

3=0 

If we apply (|7^ to V'i'^O' while simultaneously anticipating the relations f l^ . (|71|) and ((7S|l . we get the spectrum for 
Hr as 

fe 

= ^ej=eo + 2wfc = w(2fc + a + l), (77) 

coinciding with H27|) . as expected. 



5. SELF-ADJOINT EXTENSION 



In this section we investigate the self-adjoint extensions PJ|] of the radial part of the Hamiltonian. Physically this 
means that we would find all possible boundary conditions for which the radial Hamiltonian is self-adjoint. This 
would be done by an appropriate analysis of the differential operator for the radial Hamiltonian. We shall see that 
the system admits a one-parameter family of self-adjoint extensions for certain values of the system parameters. We 
start with a rederivation of the energy eigenvalues discussed in Section 2 and this method would be carried over when 
we discuss self-adjoint extensions. 

The Schrodinger's equation obeyed by the Hamiltonian of the system is given by (|47|) . The operator involving 
the radial coordinate r in the l.h.s. of H47|l satisfies the eigenvalue equation given by 



ld_ 

r dr 



I 2 2 I M " 
-— +L0 r -\ 



C 



u{r) = 2Eu{r), 



where C = {n + vY as given in l|18|) . We shall assume that ti > 0, > 0. 
In order to proceed, we make the following transformations: 



(78) 



u{r) 



t 



r e 



"X(0; 



(79) 
(80) 



and assume that /i + C > 0. In these new variables, the radial equation H78I) becomes 

X-0, 



a+l E\ 



whose solution is given by 0| 



x(r)=Af( a + l, ujT^ 



(82) 



where M denotes the confluent hypergeometric function. The above solution in general consists of an infinite series. 
However, in order for the solution to be square integrable, this series must terminate, which happens when 



a + l E _ 

where A: is a positive integer. From this, and using the definitions of b and C, we obtain that 



En,k — ^ 



(83) 



(84) 
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which is the same as (|27l) . Also, substituting (|83() in (|82|l and using the relation between confluent hypergeometric 
and Laguerre functions |13], we see that 



where the symbol (p)„ means 



- P(P + + 2) {p + n-l), {p)o = 1. 



(85) 



(86) 



The full solution of the radial equation is obtained by substituting (|85(l in (|81(l . 
We now consider the self-adjoint extensions of the operator where 



Or 



9^ 



ld_ 

r dr 



H + C 



(87) 



To begin with, we briefly recall the salient points of von Neumann's theory of self-adjoint extensions \14 

Let T be an imbounded differential operator acting on a Hilbert space H and let D{T) be the domain of T. The 
inner product of two elements, a, f3 E H is denoted by {a, (3). Let D{T*) be the set oi (p G H for which there is a 
unique G Ti with (T^,(/?) = (C,??) V ^ G D{T). For each such e D{T*), we deflne T*ip = r/. The operator T* 
then deflnes the adjoint of the operator T and D{T*) is the corresponding domain of the adjoint. The operator T is 
called symmetric or Hcrmitian iff {T(f,rj) = {ip^Trf) V (p,ri G D{T). The operator T is called self-adjoint iff T = T* 
andD{T) = D{T*). 

We now state the criterion to determine if a symmetric operator T is self-adjoint. For this purpose, let us define 
the deficiency subspaces K± = Ker(z =F T*) and the deficiency indices n±{T) = dim[ii'±]. Then T falls in one of the 
following categories: 

1) T is (essentially) self-adjoint iff (n+,n_) = (0,0). 

2) T has self-adjoint extensions iff ri+ — n^. There is a one-to-one correspondence between the self-adjoint extensions 
of T and the unitary maps from into K- . 

3) If n+ 7^ then T has no self-adjoint extensions. 

We now return to the discussion of the operator Or- This is an unbounded differential operator defined in The 
operator Or is a symmetric operator on the domain D{Or) = {¥^(0) = <(5'(0) = 0, ^p, ip' absolutely continuous, p> G 
L^(r(ir)}. Next we would like to determine if Or is self-adjoint. We shall focus on the case where a > 0. 

The deficiency indices n± are determined by the number of square-integrable solutions of the equations 

0.>±(r) = ±m±(r), (88) 

respectively, where O* is the adjoint of Or and the functions u±{r) span the deficiency subspaces K±^ respectively. 
Note that O* is given by the same differential operator as Or- From dimensional considerations we see that the r.h.s. 
of Eq. H88(l should be multiplied by a constant with dimension of length"^. We shall henceforth choose the magnitude 
of this constant to be unity by an appropriate choice of units. 
Equation (|88|l can be written as 



t) 



dt 



2uj J 



(89) 



where 



ii±(r) -r"e- — x±(r). (90) 
The solutions u±{r) of (|88|l must be square integrable. We are thus led to choose the solutions given by 



X±{r) = U {g±, a + 1, ur"^) , 



(91) 
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Figure 3. The n = and n = 1 bands, representing the parameter ranges where self-adjoint extension is possible, are drawn. 
Infinite number of such bands exists in the right half of the — v plane. 



where 



with g± = ^ T 2^^' «± 



U {g±, a + 1, Lur'^) = A 
d 



M(g±, a- 



1, ujr"^) a M [d±, 1-a, ujr'^) 

^^"^ ' T(g±)T(l-a) 



=F T- and A 

2 ' luj 



r(d±)r(a + i) ^ > r(g±)r(i-a) 

" sin(7r(a+i)) ' Consequently we obtain that 



u±(r) = r"e U [g±, a + 1, wr^) 



(92) 



(93) 



The functions u±{r) are general solutions of (|88|l with the property that as r ^ oo, u±{r) sufhciently fast such 
that they are square integrable at infinity. In order to investigate the square integrability of the functions u± (r) near 
r — 0, first note that as r ^ 0, M{g±, 6+1, ajr"^) 1. Using these, we see that as r ^ 0, 



u±\ r. 



dr. 



(94) 



where Ai,A2 and A3 are constants independent of r. From H94|l it is clear that the functions u±{r) are not square 
integrable near r = when a > 1. Thus when a > 1, we have the deficiency indices n_|_ = ?i„ — and the operator Or 
is essentially self-adjoint in the domain D[Or)- However, when < a < 1, the functions u±{r) are square integrable 
near r — Q and hence for the whole range of r. Thus, when < a < 1, we have the deficiency indices n+ = n_ = 1. 
In this case, according to von Neumann's theory, the operator Or is not self-adjoint in the domain D{Or) but admits 
a one-parameter family of self-adjoint extensions which are labeled by e*^, where z ^ R (mod 2-k). The domain of 
self-adjointness is given by D^iOr) — D{Or) ® {a{u+{r) + e*^u_(r))}, where a is an arbitrary complex number. 

Before proceeding, let us discuss the nature of the parameter ranges for which the system is either essentially 
self-adjoint or admits self-adjoint extensions. From the condition on a, we see that the self-adjoint extension exists 
when 



{)< lJL+{n + vf < 1. 



(95) 
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The eigenvalues of the angular equation are characterized by the integer n > 0. For each fixed value of n, H95|) gives a 
region on the right half of the ^ — v plane for which self-adjoint extensions exist. Let us call this region on the ^ — v 
plane as the n-th band. It is evident that such bands, corresponding to different values of n, do not overlap with each 
other. Consequently, if the values of the coupling constants /i and v are fixed on the n-th band, then the system will 
admit self-adjoint extension only if the eigenvalue of the angular equation is taken as n. It is interesting to note that 
there exists a finite gap between any two consecutive bands labeled by the integers n and n + 1. Since Eq. is 
not satisfied for any eigenvalue of the angular equation within such band gaps and also for the region ji + > 1, 
the system is essentially self-adjoint in these regions. Thus a band structure, with an infinite number of bands in the 
right half of the fi — v plane, represents the parameter ranges where the system admits self-adjoint extensions. In 
Fig. 3 we have drawn the first two of such infinite number of bands (i.e., n — and n ~ 1 band). 

In order to discuss the solutions of the eigenvalue problem, first note that in the parameter range where the system 
is essentially self-adjoint, the spectrum has already been found before, in H82I) and (|84|l . We now proceed to find the 
spectrum of Or in the domain Dz{Or) for the parameter range where the system admits self-adjoint extensions. To 
that end, first note that the solution of the eigenvalue equation (|78l) which is square integrable at infinity is given by 



u{r) = Br"e {g, a + 1, ojr^) , 

where g — — ^ and _B is a constant. In the limit r — > 0, using (|92|l . we get that 



u{r) — > AB 



r{d)T{a + l) r(5)r(i-a) 



where d = — Again, as r ^ 0, we see that 



V{a+l)\T{d+) T{d-)) r(l~a) Vr(.g+) T{g^ 
If u{r) E Dz{Or), then by comparing the coefficients of different powers of r in H97|) and H98|l we get 



fiE) 



(96) 



(97) 



(98) 



(99) 



where F + ^) = Pie 



and F {^-^ + = /92e"^^. For a given choice of the system parameters, (|^ gives 
the energy eigenvalue i? as a function of the self-adjoint parameter z. For a fixed set of system parameters, different 
choices of z lead to inequivalent quantization and to the spectrum for this model in the parameter range where it 
admits self-adjoint extension. In general, the energy E cannot be calculated analytically and has to be obtained 
numerically by plotting H99|l . a sample of which is given in Fig. 4. 

It is interesting to note that the usual A^-body Calogero model with the confining interaction leads to a similar 
radial Hamiltonian, whose self-adjoint extension has been studied before An alternate treatment of the radial 
problem in the SUSYQM framework, in terms of the self-adjoint extension of the supercharges (first order operators), 
has been given in |2^. Although the radial operators in these works share the same formal structure, the associated 
physical interpretations vary, e.g. the appearance of the band structure in the parameter space of the present model 
(as shown in Fig. 3) is absent in the case of the usual Calogero model with confining interaction. 
We conclude this section with the following observations : 

1. If r ("3^ — ^) = oo, we obtain that Ek — uo{2k + a + 1) where /c is a positive integer. For this to happen, the 
self-adjoint extension parameter z must take the value tt -|- 2(Ti and for this choice of z, we recover the usual 
eigenvalues of this system. 

2. For any other choice of z, the spectrum of Or must be obtained numerically and it is seen that the corresponding 
spectrum is not equispaced in the quantum number k. Moreover, generically there is a single negative energy 
solution. This is due to the fact that for the parameter range where this model admits self-adjoint extension, 
and for a generic value of z, SU{\^ 1) can no longer be implemented as a spectrum generating algebra. 
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Figure 4. A plot of Eq. iinnj 

using Mathematica with u) — 0.25, a = 0.25 and z = —1.5. The horizontal straight line 
corresponds to the value of the r.h.s of Eq. J99t . 



3. Although we have assumed that /i + C > 0, in general this restriction is not necessary. If this is equal to zero, the 
system can still be shown to admit a one-parameter family of self-adjoint extensions by using an analysis similar 
to the one presented here. Also, if + C is strongly negative, then this problem might require renormalization 
[i^l, which would be studied elsewhere. 

6. GENERALIZATION TO THE -BODY CASE 

Following the same procedure as for the two-body case, let us now investigate some important features of the N 
-body case. In particular, we are interested in the general structure of the spectrum which basically has the same 
structure as that shown in Fig. 2, except that for the N -body case, the Bargmann vacua depend on A'^— 1 quantum 
numbers. The problem under consideration is described by the Hamiltonian 



with A = 2v{v — 1). In Ref. |23| some other deformations of the basic inverse-square interparticle-distance- model are 
considered, but these models do not possess the SU{\, 1) symmetry. This makes them unsuitable for the treatment 
within the scope of the techniques based on the transition to the Bargmann representation. Owing to the underlying 
conformal structure of the model HIOOII . the spectrum and the eigenstates can be deduced by employing the same 
techniques as before. These are based on the fact that the Hamiltonian 1] 100(1 has some special favourable properties 
which allow us to link it with the set of decoupled oscillators that is much easier to handle with. In order to do this, 
we introduce the set of operators 



rp _ 1 Y^iV 

+ - 2 Z^i=i 

2 1^1=1 'd^ 2 ^i<i {xi-Xj)'-' 2^" UOlj 



J-0 — 2 Z^i=l dxi 



N ^ a , N 

4 



satisfying (|SJ. The connection of the Hamiltonian (|100() with the set of decoupled oscillators [29j, described by the 
operator Tq, is provided through the transformation jT)), where now T_|_, and T_, are from H101() . As a matter 
of remark, exactly this is what we mean when we are referring to the transition to the Bargmann representation. 
After this transition we are left with the equation r„A = 0, which has to be solved in order to find all Bargmann 
vacua. Of course, they will now depend on A^ — 1 quantum numbers. The remaining quantum number required for 
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the complete determination of eigenfunctions (eigenfunctions of (|100|l are fully specified with A'' quantum numbers) 
will describe excitations within each tower of states built up over the corresponding ground state (see 2). These 
excitations are governed by the collective relative radial motion of particles and are a universal feature |l5l Il7| of all 
models possessing the underlying SU (1,1) symmetry. 

In finding the solution to the equation T_A = 0, one can make a factorization 



A^Ylix^-Xj^cj), (102) 



to obtain an equation for 



One solution to this equation is of the form 




(104) 



with 



/3= • (105) 

It is clear that for each ordering of particles the solution (|104ll does not have nodes and hence it is the solution for 
the ground state. As far as we are interested in the energy and the wave function of the ground state and all vertical 
(see Fig. 1) excitations over it, the solution H104|) is all we need to deduce the spectrum and wavefunctions over the 
lowest lying vacuum in Fig. 1. In order to obtain the ground state of the model Hamiltonian (|100|l and to construct 
all its wave functions V'o.fc,o,...,o belonging to the first tower of Fig. 1, we have to use the transformation {Tj) to return 
back from the Bargmann representation describing the set of decoupled oscillators, 




V'o,M,...,o'-e--^+e"^^- ((r+)^-Aj ^ ^2 ] ]"[ _ ^^.)-e"-^-L^(2c.r+), (106) 

where L'^{2ujT+) are the associated Laguerre polynomials with 



"=2 



]-\l{N -2 + - vNf + (107) 



The corresponding part of the spectrum follows from the equation TqA = ^IS., with Tg and A having the form 
of H101|l and H102(l . respectively, and is given by 

^^o,fc,o,...,o = w(eo.o,...,o + 2k)^Lo(^ + [N - 2 + vN'^ - vNf + 4/i + 2fc^ . (108) 

Here eo is an abbreviation for eo,o,...,o designating the energy of the ground state, that is the energy of the lowest 
of the Bargmann vacua which were seen to depend on iV — 1 quantum numbers. 

Of course, there still remains the problem of finding all solutions to Eq. f ll03l) . Obviously, there are many solutions 
to Eq. ( I103|l and each of them would define one particular Bargmann vacuum. Once we find all these solutions, 
we have completed the task of integrating the model Hamiltonian (|10U|I because all excited states built over some 
particular solution A„j^,„3_...^„j^ to Eq. { I103|) are described by an associated Laguerre polynomials in the radial 
variable 2cjr-|_. So, the associated Laguerre polynomial will appear to describe excitations in each tower of states 
shown in Fig. 1 and, as we have seen, this is a common feature for all conformally invariant models. 

If we look at theproblem of integration of Eq. ( I103() more closely, we are naturally led to the question on the 
superintegrability 0| of the Hamiltonian (|100|l . This question is closely related js^ to the problem of separability 
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and even multi-separability of the corresponding Schrodinger equation. Namely, a superintegrable system is one 
that admits more integrals of motion than it has degrees of freedom, and if it is characterized by a complete set of 
commuting quadratic integrals of motion, then it is also multi-separable. This means that its Schrodinger equation 
(i.e. the Hamilton-Jacobi equation in the classical case) allows the separation of variables in more than one orthogonal 
system of coordinates. Since it is known that the N -body Calogero model [s^ is superintegrable, it would be of 
interest to see whether the extra term in HlOOfl changes anything in this respect. To get the answer to this question, 
it is most convenient to consider the three-body simplification. If we write the determining condition for Bargmann 
vacua, T_A = 0, in spherical coordinates and separate the radial part from the angular one, we are confronted with 
the set of relations 

^u+^-^u^^u^O, (109) 

A,G(M) - C)F„.„. - -i^A („„Wh.) _ = 0, (110) 

where C is the integration constant and 

G{e, dp) = 3 ^ ^ + i 2 + ^ 2 • (111) 

sin 6'(1 - sin2(/)) (sin0cos0 - cos0) (sin sin - cos 6*) 

In the above procedure it is understood that the Bargmann vacua depend on 2 quantum numbers, as is readily 
expected from the general consideration, and is separated as 

An,m = u(r)F„,„(0,<^). (112) 



While Eq. ( 1109(1 is easily integrated to u ~ ^ ^}jg angular equation remains a much more difficult 

problem to solve. As has already been said, this problem amounts to the question of superintegrability of the 3 - 
body version of the Hamiltonian (|100|l and consequently to the problem of separability of Eq. H110|) . Nevertheless, 
we can deduce the whole spectrum by relying only on the expression 1)112(1 and on the fact that the spectrum has to 
become a Calogero one in the limit when the parameter /i approaches zero. With these observations in mind, the 
whole spectrum for the 3-body variant of the Hamiltonian ((100() is readily found to be 

En,k,rn ^ UJ + 2k + ^^1+4(3 V + n + 3m)(l + 3t^ + ?i + 3m) + . (113) 



7. CONCLUSION 



We have considered a conformally invariant deformation of the quantum Calogero model with the special emphasis 
on the deformation of the two-body model. Owing to the fact that the model under consideration possesses SU{1, 1) 
symmetry, we were able to apply the techniques based on the correspondence between some particular conformally 
invariant model and the set of decoupled oscillators, the procedure which is referred to as the Bargmann representation 
analysis. This has provided us with the possibility to construct the creation and annihilation operators acting on the 
Fock space, thereby allowing us to investigate the dynamical structure of the problem. 

We have also analyzed the self-adjoint extensions of the radial Hamiltonian of the two-body problem and have found 
the region in the parameter space where the system admits a one-parameter family of self-adjoint extensions. In the 
situations where the system admits self-adjoint extensions, SU(1, 1) can no longer in general be implemented as the 
spectrum generating algebra and the corresponding spectrum is not equispaced in the quantum number k. However, 
for a special value of the self-adjoint extension parameter, SU{1, 1) can be recovered as the spectrum generating 
algebra, in which case the spectrum becomes equispaced in fc. It is plausible that the angular part of the Hamiltonian 
given in ((49(1 might also admit self-adjoint extensions, anaologous to the angular Hamiltonian of the 3-body Calogero 
model as discussed in Isil. 
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At the end, we have carried out the generahzation to include the N— body problem as well. Here we have found 
the ground state, which is the lowest of all Bargmann vacua, and excitations over it, together with the corresponding 
spectrum. Also, we have found that all excitations are basically the same, all of them having the origin in the collective 
relative radial motion of the particles and are described by the associated Laguerre polynomials in the collective radial 
variable. It appears that this is, in fact, a common feature of all models possessing the underlying SU{1, 1) symmetry. 
The only problem, but in no case the simple one, that has been left is to find all Bargmann vacua, that is to find 
all solutions to Eq. Ij 103(1 . We have seen that this problem is closely related to the problem of the superintegrability 
of the Hamiltonian in question. To somehow clear up the situation, we have made the simplification to consider the 
three-body problem. Although we have not found how all Bargmann vacua look like in the three-body problem, we 
were able to deduce the complete spectrum for this case. All that was needed was the information on the underlying 
conformal invariance of the model and the explicit form of the radial part of the Bargmann vacua wave functions, 
together with the note that the Calogero spectrum has to be reproduced in a smooth limit when the deformation 
parameter tends to zero. We hope that the issues that are left unresolved, as is the case with the structure of the 
Bargmann vacua in a general N— body problem, will be addressed in the near future, at least for the three-body 
case. 
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